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Abstract:  This  paper  presents  a  detailed  investigation  of  the  motion  of 
individual  micro-particles  in  a  moderately- viscous  liquid  in  direct  response 
to  a  local,  laser-induced  temperature  gradient.  By  measuring  particle 
trajectories  in  3D,  and  comparing  them  to  a  simulated  temperature  profile,  it 
is  confirmed  that  the  thermally-induced  particle  motion  is  the  direct  result  of 
thermophoresis.  The  elevated  viscosity  of  the  liquid  provides  for  substantial 
differences  in  the  behavior  predicted  by  various  models  of  thermophoresis, 
which  in  turn  allows  measured  data  to  be  most  appropriately  matched  to  a 
model  proposed  by  Brenner.  This  model  is  then  used  to  predict  the  effective 
force  resulting  from  thermophoresis  in  an  optical  trap.  Based  on  these 
results,  we  predict  when  thermophoresis  will  strongly  inhibit  the  ability  of 
radiation  pressure  to  trap  nano-scale  particles.  The  model  also  predicts  that 
the  thermophoretic  force  scales  linearly  with  the  viscosity  of  the  liquid,  such 
that  choice  of  liquid  plays  a  key  role  in  the  relative  strength  of  the 
thermophoretic  and  radiation  forces. 
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1.  Introduction 

It  has  long  been  known  that  particles  immersed  in  a  gas  or  liquid  will  migrate  in  response  to  a 
thermal  gradient,  an  effect  known  as  thermophoresis.  In  liquids  this  leads  to  phenomena  such 
as  the  Soret  effect  [1,2],  in  which  mixtures  tend  to  separate  in  response  to  a  thermal  gradient. 
In  gases  the  origins  of  thermophoresis  are  reasonably  well -understood,  and  the  relationship 
between  thermal  gradient  and  induced  motion  can  be  related  directly  to  the  kinetic  theory  of 
gases  [3].  For  liquids,  however,  there  are  numerous  contradictory  theoretical  descriptions  [3- 
11],  partially  due  to  the  lack  of  experimental  data  on  single -particle  motion  [3,7],  and  partially 
to  the  fact  that  thermo -diffusion,  having  to  do  with  the  concentration  of  particles  in  a  liquid,  is 
most  easily  treated  empirically  [1]. 

In  the  context  of  laser  trapping  and  manipulation  of  small  particles,  the  thermal  force 
induced  by  optical  absorption  within  a  particle ,  and  the  resulting  temperature  gradient  (known 
as  the  indirect  photophoretic  force),  has  been  well  explored  for  particles  suspended  in  a  gas 
[12].  Indeed,  this  force  forms  the  basis  for  the  “radiometric  barrier”  which  must  be  overcome 
to  trap  particles  at  low  gas  pressure  [13].  In  liquids,  laser-induced  thermo -diffusion  of  large 
collections  of  nanoparticles  has  also  been  studied,  but  primarily  in  the  context  of  thermally - 
induced  variations  in  steady-state  particle  concentration  (i.e.  the  Soret  effect)  rather  than  the 
motion  of  individual  particles.  Empirical  and  theoretical  models  describing  thermo -diffusion, 
however,  tend  to  differ  from  the  observed  motion  of  individual  particles  in  a  liquid  [3,14]. 
Studies  of  laser-induced  thermophoresis  of  individual  particles  in  liquids  have  also  been 
performed  [10,15,16],  but  quantitative  results  are  generally  limited  to  aqueous  solutions,  and 
thus  provide  limited  data  for  the  evaluation  of  theory.  Such  works  also  rely  on  very  weakly- 
focused  optical  beams,  quite  different  from  that  of  typical  optical  tweezers. 
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For  traditional  (radiation-pressure-based)  optical  traps,  the  thermophoresis  of  individual 
particles  in  liquids  has  received  little  attention.  This  is  due  at  least  in  part  to  the  dominance  of 
radiation  pressure  over  thermal  forces  in  typical  optical  trapping  experiments,  such  that 
thermophoretic  migration  is  not  easily  observed  in  a  trap.  A  lack  of  theoretical  understanding 
of  thermophoresis  in  liquids,  however,  has  also  undoubtedly  played  a  part.  A  consequence  is 
that  standard  theoretical  models  for  forces  in  an  optical  trap  do  not  account  for  a 
thermophoretic  force  [17,18].  Given  that  thermophoresis  is  a  proven  method  for  manipulating 
groups  of  nano-scale  particles  in  liquids,  whereas  radiation  pressure  struggles  on  the  nano¬ 
scale  [10,19],  it  is  important  to  bridge  the  gap  between  the  two  phenomena. 

This  paper  attempts  to  shed  some  light  on  the  subject,  by  characterizing  the  laser-induced 
thermophoresis  of  individual  microspheres  in  a  moderately-viscous  liquid  (sixty  times  more 
viscous  than  water).  The  elevated  viscosity  of  the  liquid  leads  to  substantial  differences  in  the 
behavior  predicted  by  various  models  of  thermophoresis,  thus  allowing  measured  data  to  be 
most  appropriately  matched  to  a  theoretical  model  proposed  by  Brenner.  This  model  is  then 
used  to  predict  the  effective  force  resulting  from  thermophoresis  in  an  optical  trap.  Based  on 
these  results,  we  predict  when  thermophoresis  will  strongly  inhibit  the  ability  of  radiation 
pressure  to  trap  nano-scale  particles.  The  model  also  predicts  that  the  thermophoretic  force 
scales  linearly  with  the  viscosity  of  the  liquid,  such  that  choice  of  liquid  plays  a  key  role  in  the 
relative  strength  of  the  thermophoretic  and  radiation  forces. 

2.  Experiment 
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Fig.  1 .  Schematic  diagrams  of  (a)  the  optical  system  used  to  manipulate  and  track  particles  in  a 
liquid,  and  (b)  the  liquid-filled,  square  cross-section,  glass  capillary  sample  cell.  The  cell  was 
held  suspended  in  air  by  a  rigid  clamp,  located  2.0  cm  above  the  laser  focus  in  the  positive  y- 
direction.  Gravity  was  directed  along  the  negative  y-direction. 


A  schematic  diagram  of  the  experimental  setup  used  to  study  laser  thermophoresis  in  liquids 
is  shown  in  Fig.  1 .  The  system  utilized  a  simple  single-beam  optical  trap,  formed  by  launching 
the  output  of  an  808  nm  distributed  feedback  laser  diode  through  a  collimating  lens,  reflecting 
it  off  of  a  dichroic  mirror,  and  focusing  to  a  spot  with  an  infinity -corrected,  0.66  numerical 
aperture  objective  lens.  At  the  focus  of  the  trap  was  placed  a  sample  cell,  consisting  of  a 
liquid-suspension-filled,  borosilicate  glass  capillary  with  square  cross-section  (50  pm  inner 
width,  25  pm  wall  thickness)  and  5  cm  length.  The  cell  was  oriented  such  that  the  capillary 
wall  was  normal  to  the  incident  beam,  and  its  long  axis  was  vertical  (parallel  to  the  force  of 
gravity).  It  was  held  suspended  in  air  by  a  rigid  clamp,  which  first  contacted  the  cell  2  cm 
above  the  beam  focus.  The  cell  extended  roughly  1  cm  below  the  beam  focus.  Due  to  the  flat 
walls  of  the  capillary,  and  the  narrow  thickness  of  fluid  within,  the  position  of  individual 
particles  in  the  liquid  could  be  easily  observed  along  two  orthogonal  directions.  This  allowed 
two  CCD  cameras,  shown  in  Fig.  1,  to  image  and  track  particle  positions  in  3D,  one 
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measuring  particle  location  (x,y),  and  the  other  (y,z).  For  the  (y,z)  view,  a  lower,  0.28 
numerical  aperture  lens  was  used  to  provide  greater  depth  of  field.  3D  video  was  acquired  by 
triggering  the  exposure  of  both  cameras  simultaneously,  at  30  frames  per  second  with  a 
common  signal,  and  recording  the  images  on  a  computer.  Illumination  of  the  cell  was 
performed  using  a  pair  of  white  light  emitting  diodes  (LEDs),  which  were  chosen  to  prevent 
sample  heating  due  to  infrared  absorption  of  the  illumination  source.  The  divergence  of  the 
trap  beam  was  measured  to  be  36°  in  air,  by  profiling  the  beam  intensity  with  a  z-scanned 
mirror.  The  optical  power  of  the  trap  beam  was  also  measured  to  be  38  mW  in  air.  After 
accounting  for  Fresnel  reflections  in  the  sample  cell  the  power  was  estimated  to  be  35  mW  in 
the  liquid. 

In  order  to  investigate  particle  motion  due  to  laser-induced  thermophoresis,  it  was 
necessary  to  limit  the  relative  strength  of  the  radiation  pressure  produced  on  the  particles  in 
the  sample  cell.  One  option  would  have  been  to  use  a  highly-absorptive  liquid  (or  wavelength) 
to  enhance  the  temperature  gradient  relative  to  the  radiation  pressure  on  the  particle.  An 
alternative  was  to  reduce  the  radiation  pressure  on  the  particle,  without  significantly  altering 
the  optical  absorption  of  the  liquid.  The  latter  was  chosen  for  this  experiment,  because  it  was 
believed  to  provide  a  better  representation  of  the  thermal  conditions  in  a  typical  optical 
trapping/manipulation  experiment. 

In  order  to  reduce  the  radiation  pressure  on  the  particles,  without  resorting  to  extremely 
small  particles  which  would  be  more  susceptible  to  Brownian  motion,  it  was  necessary  to 
reduce  the  refractive  index  contrast  An  =  np-nh  between  the  particle  and  liquid,  where  np  and 
rii  are  the  refractive  indices  of  the  particle  and  liquid,  respectively.  This  was  accomplished  by 
mixing  dry  silica  microspheres  (np  =  1.45  [20])  with  99.5%  pure  propylene  glycol  (1,2 
propanediol)  ( =  1.43  [20]),  forming  a  low-concentration,  low  index  contrast  (An  =  0.02) 
suspension.  The  suspension  was  then  placed  in  the  sample  cell  via  capillary  action.  In  the 
small-particle  limit,  radiation  pressure  can  be  separated  into  gradient  and  scattering  force 
components  which  scale  roughly  as  An  and  (An)2,  respectively  [17].  Thus,  compared  to 
polystyrene  spheres  in  water  (An  =  0.25)  [21],  the  radiation  pressure  in  the  experiment  should 
have  been  substantially  reduced.  Furthermore,  since  the  optical  absorption  of  propylene  glycol 
is  similar  to  that  of  water  [22],  as  is  its  thermal  diffusivity,  the  temperature  distribution  within 
the  cell  should  have  been  comparable  to  that  of  a  water-filled  cell.  An  additional  advantage  of 
propylene  glycol  was  that  it  provided  a  relatively  high  viscosity,  0.055  kg/m-s,  which  served 
to  dampen  Brownian  motion  of  the  particles  in  the  experiment,  and  to  suppress  convective 
flow. 

Figure  2  presents  a  time-lapsed  video  (Media  1)  of  the  laser-induced  motion  of  two  silica 
microspheres  in  propylene  glycol.  Time  in  the  video  has  been  scaled  by  a  factor  of  30,  such 
that  one  second  in  the  video  corresponds  to  thirty  seconds  in  the  experiment.  Also  shown  are 
curves  representing  the  1/e2  intensity  limits  of  the  optical  beam,  as  calculated  by 
approximating  a  Gaussian  beam  profile  of  the  form 


(1) 


where 


(2) 
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Fig.  2.  Time-lapsed  video  (Media  1)  of  the  motion  of  two  silica  microspheres,  (A)  and  (B), 
suspended  in  propylene  glycol,  in  response  to  a  focused  laser  beam.  Each  second  in  the  video 
corresponds  to  thirty  seconds  in  the  experiment.  The  1/e2  intensity  limits  of  the  laser  beam, 
calculated  using  Eqs.  (1)  and  (2),  have  been  overlaid  in  gray.  The  thicknesses  of  the  cell  walls 
and  liquid  region  between  were  25  pm  and  50  pm,  respectively.  Microsphere  (A)  falls  into  the 
beam,  becomes  trapped  in  (x,y),  and  is  subsequently  pushed  along  the  beam  24  pm  before 
being  released  from  the  2D  trap.  Microsphere  (B)  falls  past  the  beam  without  becoming 
trapped,  but  experiences  a  similar  push  when  closest  to  the  beam  axis.  Particle  size  is 
exaggerated  due  to  imperfect  focus. 


_  jo  w  \ 

7inl  tan  60l  n  tan  6{)a 


(3) 


P  is  the  optical  power,  X0  is  the  vacuum  wavelength,  co0  is  size  of  the  beam  waist,  and  00i  and 
60a  are  the  beam  divergences  in  liquid  and  air,  respectively  [23].  Using  the  measured  beam 
divergence  in  air,  the  waist  size  was  estimated  to  be  0.4  pm.  The  location  of  the  beam  waist  in 
the  video,  which  corresponded  to  the  origin  of  the  coordinate  system,  was  calibrated  by 
moving  the  cell  in  z  with  a  precision  micrometer  to  focus  on  different  microspheres  with 
camera  1,  and  noting  their  locations  in  the  camera  2  view.  Since  the  beam  waist  coincided 
with  the  image  focus  (in  air),  this  calibration  was  taken  as  a  reasonable  estimate  of  the  beam 
waist  position.  The  location  of  the  beam  axis  (v,y)  =  (0,  0)  was  readily  determined  by  imaging 
the  reflected  beam  profile  with  camera  1.  With  this  calibration,  the  central  axis  of  the  capillary 
tube  coincided  with  the  coordinates  ( x,z )  =  (-1.6,14)  pm. 

Figure  2  demonstrates  that  the  microspheres  initially  fell  in  the  negative  y  direction,  which 
was  due  to  a  combination  of  gravitational  force  (resisted  by  viscous  drag)  and  a  creeping 
laminar  flow  resulting  from  evaporation  from  the  open  ends  of  the  capillary.  Upon  nearing  the 
laser  beam,  microsphere  A  was  pulled  rapidly  into  the  center  of  the  beam,  becoming  trapped 
in  x  and  y,  and  subsequently  pushed  along  the  beam  a  distance  of  24  pm  before  being  released 
from  the  trap  and  resuming  its  downward  motion.  In  the  same  timeframe,  microsphere  B 
passed  by  the  optical  beam,  with  no  noticeable  change  in  direction  except  when  nearest  the 
optical  beam. 

Figure  3  plots  the  positions  of  microspheres  A  and  B  in  3D  as  a  function  of  time,  which 
were  determined  by  particle  tracking  analysis  of  the  image  sequences  from  both  cameras.  As 
shown,  microsphere  A  was  initially  pulled  into  the  beam  and  pushed  rapidly  towards  the 
waist,  and  then  upon  passing  the  waist  exhibited  slow  but  relatively  constant  velocity  along 
the  direction  of  the  beam.  Past  the  beam  waist,  it  moved  with  an  average  velocity  0.12  pm/s 
over  a  19  pm  distance,  while  at  the  same  time  remaining  trapped  in  v  and  y.  Microsphere  B, 
on  the  other  hand,  which  was  offset  from  the  beam  axis  by  -5.5  pm  in  the  x-direction,  and 
therefore  never  trapped,  exhibited  a  similar  0.11  pm/s  velocity  when  within  8  pm  of  the  beam 
axis.  This  was  despite  the  fact  that  microspheres  A  and  B  never  came  in  contact  with  each 


#144014  -  $15.00  USD  Received  14  Mar  2011;  revised  18  Apr  2011;  accepted  1  May  2011;  published  13  May  2011 
(C)  201 1  OSA  23  May  201 1  /  Vol.  19,  No.  1 1  /  OPTICS  EXPRESS  10575 


Fig.  3.  Positions  of  (a)  microsphere  (A)  and  (b)  microsphere  (B),  as  a  function  of  time.  The 
origin  of  the  coordinate  system  coincided  with  both  the  beam  axis  and  the  beam  waist. 
Microsphere  (A)  became  trapped  in  2D,  at  (x,y,z)  =  (0,0,z),  and  was  then  pushed  in  the  positive 
z-direction.  Microsphere  (B)  was  not  trapped,  but  was  also  pushed  in  the  positive  z-direction 
when  within  8  pm  of  the  beam  axis. 


other,  as  their  radii  were  1.6  jam  and  1.7  jam,  respectively  (as  measured  using  camera  1),  but 
their  closest  center-to-center  spacing  was  7.5  pm.  The  difference  in  y- velocity  between  the 
two  microspheres,  despite  their  similar  size,  can  be  attributed  to  the  creeping  fluid  flow  in  the 
negative  y-direction,  which  would  be  expected  to  exhibit  a  parabolic  velocity  profile.  This 
flow  was  estimated  to  have  a  peak  y-velocity  of  0.2  pm/s,  by  tracking  the  y-velocities  of 
various  microspheres  of  different  size  and  location  within  the  cell,  and  correcting  for  terminal 
velocity  due  to  gravity. 

3.  Analysis 

The  resulting  z-velocity  of  microsphere  A  is  plotted  versus  position  along  the  beam  axis  in 
Fig.  4(a).  This  curve  was  calculated  by  first  fitting  the  position  versus  time  data  to  a  smoothed 
spline  and  then  differentiating,  which  helped  to  reduce  the  effects  of  noise  (though  noise  due 
to  mechanical  vibrations  and  particle  tracking  error  was  accentuated  by  differentiation,  the 
overall  trend  is  clear  in  both  Figs.  3  and  4).  As  shown,  microsphere  A  experienced  a  rapid 
acceleration  towards  the  beam  waist  (at  z  =  0),  which  then  all  but  vanished  upon  reaching  the 
waist.  Such  behavior  is  qualitatively  consistent  with  a  3D  optical  trap,  which  pushes  a  particle 
towards  an  equilibrium  point  located  just  beyond  the  beam  waist  via  radiation  pressure  [17]. 
The  motion  observed  beyond  and  away  from  the  waist,  in  the  region  z>0 ,  was,  however, 
inconsistent  with  a  stable  3D  trap. 


Fig.  4.  Measured  z-velocity  versus  axial  position,  compared  to  (a)  the  z-velocity  predicted  by 
radiation  pressure,  and  (b)  the  local  intensity  of  a  Gaussian  beam.  In  (a),  the  curve  labeled 
“ideal”  is  based  on  rigorous  calculations  for  an  ideal  beam,  while  the  curve  labeled  “aberrated” 
is  a  qualitative  estimate  of  the  impact  of  weak  spherical  aberration. 
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3.1  Radiation  pressure 

In  order  to  both  illustrate  that  the  observed  motion  was  inconsistent  with  a  stable  3D  trap  and 
to  aid  in  further  analysis,  a  series  of  theoretical  calculations  of  the  optical  force  due  to 
radiation  pressure  were  performed  for  microsphere  A,  using  an  approach  based  on  Nieminen 
[24-26] .  Since  details  of  these  calculations  are  not  of  immediate  importance,  the  methodology 
has  been  summarized  in  Appendix  A.  Here  it  suffices  to  note  that  the  calculations  were 
performed  assuming  a  0.4  pm  beam  waist,  1.6  pm  particle  radius,  35  mW  optical  power,  and 
liquid  and  particle  refractive  indices  of  1.43  and  1.45,  respectively.  The  calculated  optical 
force  Frp  on  the  microsphere  was  subsequently  converted  to  a  velocity  vRP  using  Stokes’  drag 
relation 


v 


RP 


6nrjR 


(4) 


where  rj  is  the  dynamic  viscosity  of  the  liquid,  0.055  kg/m-s  [27],  and  R  the  radius  of  the 
sphere,  1.6  pm. 

The  velocity  predicted  for  micro  sphere  A  as  a  function  of  position  along  the  beam  axis  is 
plotted  in  Fig.  4(a).  As  shown,  in  the  region  z  <0  the  measured  velocity  was  in  reasonably 
good  agreement  with  that  predicted  due  to  radiation  pressure,  particularly  in  terms  of  its 
magnitude.  On  the  opposite  side  of  the  waist,  however,  agreement  was  extremely  poor.  One 
could  ask  whether  this  discrepancy  was  due  to  an  incorrect  assumption  of  a  perfect, 
aberration-free  optical  beam  in  the  radiation  pressure  calculations.  It  is  known  that  even  slight 
spherical  aberration  along  the  optical  path,  which  tends  to  blur  the  focused  spot,  will 
significantly  weaken  the  axial  restoring  force  on  the  positive  z  side  of  the  waist  [28,29].  In 
previous  studies  it  has  been  shown  that  for  weak  spherical  aberration,  the  forward  force  in  the 
z<  0  region  remains  largely  unchanged,  while  the  reverse  force  for  z>0  is  significantly 
reduced.  Based  on  data  presented  in  such  works,  we  have  included  a  qualitative  estimate  of 
how  the  velocity  of  the  microspheres  might  have  been  impacted  by  slight  spherical  beam 
aberration.  This  serves  to  illustrate  that  the  axial  restoring  force  for  z>0  due  to  radiation 
pressure  may  have  been  significantly  weaker  than  that  predicted  for  the  ideal  beam.  This  is  an 
important  point,  because  the  discrepancy  between  radiation  pressure  calculations  and 
experiment  in  the  region  z>0 ,  both  in  magnitude  and  direction,  force  us  to  consider  additional 
contributions  to  the  motion  of  the  particle,  which  must  ultimately  be  reconciled  with  the 
restoring  force  due  to  radiation  pressure.  Spherical  aberration  was  also  to  be  expected  in  the 
experiment,  since  the  lens  was  aberration-corrected  for  a  170  pm  thick  coverslip,  but  the 
experimental  focus  was  only  39  pm  into  the  sample  cell.  In  the  following  subsection  it  will  be 
shown  that  this  qualitative  estimate  leads  to  a  consistent  description  of  the  observed  motion,  in 
which  radiation  pressure  dominates  for  z<0  and  thermophoresis  dominates  for  z>8  pm,  despite 
reasonable  uncertainty  due  to  aberration  in  the  region  between. 

Figure  4(b)  compares  the  measured  z-velocity  of  microsphere  A  to  the  normalized 
intensity  along  the  beam  axis,  defined  as  Inorm  =  1(0,0, z)/ 1( 0,0,0).  As  shown,  the  rapid  1/z 
roll-off  in  intensity  for  large  z  was  very  dissimilar  to  the  relatively  constant  velocity  observed 
from  z  =  0  to  18  pm.  Based  on  this  simple  observation,  a  direct  relationship  between  the  local 
intensity  and  the  observed  motion  in  the  range  z  >  0  could  be  ruled  out.  This  was  important, 
because  it  allowed  a  wide  range  of  effects  to  be  confidently  ruled  out,  including  indirect 
photophoresis,  momentum  transfer  due  to  optical  absorption  within  the  particle,  and  radiation 
pressure,  all  of  which  would  have  produced  forces  proportional  to  either  the  local  intensity  or 
its  gradient.  Based  on  this  result,  it  was  necessary  to  consider  an  indirect  relationship  between 
the  local  intensity  and  the  observed  motion  in  the  region  z>  0. 

3.2  Thermophoresis 

A  reasonable  starting  point  was  to  consider  the  impact  of  optical  absorption  within  the  liquid. 
Optical  absorption  converts  optical  energy  to  heat  within  a  material,  according  to  the  relation 
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Q  =  2al 


(5) 


where  Q  is  the  change  in  heat  density  per  unit  time,  2a  is  the  optical  power  absorption 
coefficient  of  the  material,  and  I  is  the  local  intensity.  This  in  turn  is  related  to  the  local 
temperature  T,  through  the  heat  equation  [30] 

v*(*vr)+G  =  o,  (6) 

and  the  boundary  conditions 

Tx=T2  (7) 

kxh  •  V7]  =  k2n  •  VT2  (8) 

for  an  interface  between  materials  1  and  2.  Here  k  represents  the  thermal  conductivity  of  the 
material,  and  n  is  the  normal  unit  vector  at  the  interface.  Convective  heat  transfer  has  been 
neglected  in  Eq.  (6),  due  to  the  extremely  low  fluid  velocities  observed  in  the  experiment. 
Steady-state  conditions  have  also  been  assumed. 

In  order  to  solve  Eqs.  (5)  though  (8),  and  obtain  the  temperature  distribution  within  the 
sample  cell,  simulations  were  performed  using  commercially  available  3D  finite  element 
software.  Material  and  beam  properties  used  in  the  simulations  are  provided  in  Table  1  (the 
power  absorption  coefficient  for  propylene  glycol  [22]  was  consistent  with  independent 
measurements  performed  using  a  spectrophotometer)  .  The  air-glass  boundary  was  modeled 
using  the  heat  flux  boundary  condition 

kcn»VT  =  -h(T-Tw),  (9) 


which  approximated  convective  cooling  of  the  cell  by  the  surrounding  air.  The  ends  of  the 
cell,  located  at  y  =  +1.6  mm  in  the  simulation,  were  modeled  as  thermal  insulators.  For 
reasons  to  be  given  shortly,  the  temperature  profile  in  the  cell  was  modeled  without  taking  the 
presence  of  any  microspheres  into  account. 


Table  1.  Properties  Used  for  Thermal  Simulations 


Symbol 

Value 

Description 

ho 

0.808  pm 

Vacuum  wavelength 

C0o 

0.4  pm 

Beam  waist  size 

ni 

1.43 

Refractive  index  of  liquid  [20] 

P 

35  mW 

Optical  power  in  liquid 

2ai 

2.9  m"1 

Power  absorption  coefficient  of  liquid  [22] 

2ac 

0.1  m_1 

Power  absorption  coefficient  of  cell  [20] 

k 

0.196  W/m-K 

Thermal  conductivity  of  liquid  [27] 

kc 

1.0  W/m-K 

Thermal  conductivity  of  cell  [20] 

h 

40  W/nr-K 

Heat  transfer  coefficient  of  air  [31] 

Too 

293  K 

Ambient  temperature 
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Fig.  5.  Simulated  (a)  temperature  T  and  (b)  magnitude  of  the  temperature  gradient  in  the  z- 
direction  l<2TA%l,  in  the  y  =  0  plane.  The  black  square  marks  the  liquid-glass  boundary,  and  had 
sides  50  pm  long.  The  peak  temperature  was  0.60  K  above  ambient.  The  white  regions  of  (b) 
represent  regions  where  \dT/dz\  exceeds  13,000  K/m. 


Figure  5(a)  plots  the  resulting  temperature  distribution  in  the  sample  cell  across  the  xz 
plane  (y  =  0).  The  magnitude  of  the  z-component  of  the  temperature  gradient,  \dY/di\,  is  also 
plotted  in  Fig.  5(b).  Note  that  the  peak  simulated  temperature  in  the  cell  was  only  0.60  K 
above  the  ambient  temperature.  This  is  significant  because  it  indicates  that  temperature 
variations  in  the  cell  only  weakly-perturbed  the  material  properties  in  the  experiment. 

There  are  a  number  of  theoretical  models  that  relate  the  temperature  gradient  in  a  liquid  to 
the  particle  velocity  induced  by  thermophoresis  [3-11].  The  physics  behind  each  can  in 
principle  be  applied  to  an  arbitrary  temperature  profile.  However,  calculations  are  greatly 
simplified  if  one  can  approximate  a  spatially  uniform  temperature  gradient,  VT0 ,  where  T0  is 
the  temperature  in  the  absence  of  the  particle.  Under  such  a  condition,  each  model  yields  an 
expression  of  the  form 


(10) 


for  the  thermophoretic  velocity  of  a  spherical  particle  in  a  liquid.  The  thermal  diffusivity  DT 
takes  different  forms  for  the  different  models,  and  will  be  discussed  in  greater  detail  in  the 
following  section.  The  present  emphasis  is  on  comparing  the  form  of  Eq.  (10)  to  measured 
data. 

In  doing  so,  it  is  important  to  first  note  that  the  simulated  temperature  profile  in  Fig.  5(a) 
does  not  strictly  adhere  to  the  condition  implicit  in  Eq.  (10)  of  a  spatially-uniform  temperature 
gradient.  However,  for  all  models  of  thermophoresis  to  be  considered,  it  is  sufficient  that  VT0 
instead  vary  little  over  the  extent  of  the  particle  to  appropriately  apply  Eq.  (10).  This  condition 
is  reasonably  well-satisfied  for  a  3.2  pm  diameter  particle  near  the  center  of  the  cell,  and  thus 
the  latter  portion  of  the  trajectory  of  microsphere  A.  Closer  to  the  beam  waist,  however, 
variations  in  the  temperature  gradient  over  the  extent  of  this  microsphere  were  expected  to 
introduce  error  when  applying  Eq.  (10).  For  microsphere  B  the  rapid  variation  in  the  x  and  y 
components  of  the  temperature  gradient  (nonexistent  for  microsphere  A)  would  also  have 
introduced  error  when  using  Eq.  (10).  Thus,  in  the  following  it  is  understood  that  the 
application  of  Eq.  (10)  to  the  temperature  profile  of  Fig.  5(a)  provides  approximate  results, 
except  for  the  latter  portion  of  the  trajectory  of  microsphere  A  where  error  should  be 
negligible.  Furthermore,  the  presence  of  each  microsphere  in  the  liquid  would  have  altered  the 
temperature  profile,  not  only  due  to  heat  conduction,  which  is  accounted  for  in  DT ,  but  also  by 
altering  the  local  heat  density  Q ,  due  to  the  negligible  optical  absorption  of  the  particle.  This 
effect  was  most  pronounced  when  the  sphere  is  at  the  beam  focus,  but  elsewhere  relatively 
weak. 
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Fig.  6.  (a)  Comparison  of  the  measured  z-velocity  of  microsphere  (A)  along  the  beam  axis  with 
predicted  velocities  due  to  thermophoresis,  radiation  pressure,  and  the  combined  effect  of 
thermophoresis  and  radiation  pressure,  and  (b)  the  predicted  thermophoretic  z-velocity  of 
microsphere  (B)  versus  y  position.  The  radiation  pressure  curve  in  (a)  matches  the  estimated 
plot  for  an  aberrated  beam  in  Fig.  4(a). 

The  thermophoretic  velocity  of  microsphere  A  at  different  points  along  the  beam  axis  was 
calculated  using  Eq.  (10)  and  the  simulated  temperature  gradient,  and  is  plotted  in  Fig.  6(a). 
For  this  calculation  a  thermal  diffusivity  of  DT  =  22  pm2/s-K  was  assumed  in  order  to  provide 
a  best  fit  to  the  measured  data,  which  will  later  be  shown  to  be  consistent  with  theory.  Note 
that  velocity  predictions  have  been  omitted  very  close  to  the  beam  waist,  where  VT0  varied 
substantially  over  the  extent  of  a  3.2  pm  diameter  microsphere,  and  where  the  (relatively  low) 
absorption  of  the  microsphere  would  have  most  impacted  the  thermal  profile.  Figure  6(a)  also 
plots  the  measured  velocity  of  microsphere  A,  and  the  estimated  velocity  due  to  radiation 
pressure  from  Fig.  4(a)  for  an  aberrated  beam.  The  most  important  result  of  Fig.  6(a)  is  that 
thermophoresis  accounts  for  the  observed  nonzero  velocity  of  the  particle  far  from  the  beam 
waist,  where  radiation  pressure  falls  out  of  consideration.  Also  significant  is  that  the  combined 
effect  of  radiation  pressure  and  thermophoresis,  plotted  as  the  sum  of  the  two  velocities  in  the 
figure,  adequately  describes  the  relatively  constant  velocity  observed  in  the  experiment,  over  a 
range  where  the  local  intensity  varied  by  over  an  order  of  magnitude.  Although  the  radiation 
pressure  in  Fig.  6(a)  was  merely  an  estimate  (for  z>0),  the  resulting  agreement  between  theory 
and  experiment  provides  justification  for  its  use. 

Figure  6(a)  presents  compelling  evidence  that  the  particle  motion  observed  in  the  region  z 
>  3  pm  was  dominated  by  thermophoresis,  but  it  does  not  constitute  conclusive  proof.  The 
case  can  be  strengthened,  however,  by  considering  the  trajectory  of  microsphere  B.  This 
sphere  fell  continuously  along  a  path  that  remained  constant  in  x  (x  =  -5.5pm),  but  varied 
slightly  in  z  (z  =  13  to  17  pm),  due  to  a  positive  z-velocity  that  was  only  discernable  over  the 
limited  range  y  =  -6  to  6  pm.  For  the  purpose  of  comparison,  Fig.  6(b)  plots  the 
thermophoretic  z-velocity  predicted  for  a  similar  path,  along  the  line  (x,  z)  =  (-5.5,  15)  pm. 
This  figure  shows  that  thermophoresis  should  only  have  been  significant  over  a  very  similar 
range  of  y.  The  slight  discrepancy  between  the  predicted  8  pm  FWHM  interaction  region  and 
the  measured  12  pm  range  was  easily  attributable  to  the  variation  in  VT  over  the  3.4  pm 
diameter  of  the  sphere.  Such  agreement  between  the  measured  particle  velocity  and  the  local 
temperature  gradient,  along  multiple  axes,  provides  strong  support  for  the  conclusion  that  the 
observed  motion  in  the  region  z  >  3  pm  was  due  to  thermophoresis. 

3.3  Transfer  of  optical  momentum  to  the  liquid 

For  completeness,  it  is  also  worth  noting  that  optical  momentum  absorbed  by  the  liquid  can 
induce  convective  flow,  an  effect  known  to  be  significant  under  certain  (but  different) 
circumstances  [32].  Given  that  the  linear  momentum  flux  of  the  optical  beam  is  given  by 
Sn/c,  where  S  is  the  Poynting  vector  and  c  the  vacuum  speed  of  light,  it  can  easily  be  shown 
that  absorption  of  optical  momentum  by  the  liquid  produces  to  the  local  force  density 
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fabs=i2aIn,/C)z ■  C11) 

The  convection  currents  generated  by  this  force  density  within  the  sample  cell  were 
estimated  by  performing  3D  fluid  dynamics  simulations,  based  on  the  Navier-Stokes 
equations,  using  commercial  finite  element  modeling  software.  Figure  (7)  plots  the  resulting 
z-velocity  profile  of  the  liquid  in  the  y  =  0  plane.  As  shown,  the  peak  velocity  was  0.0036 
pm/s,  which  was  less  than  3%  that  measured  in  the  experiment.  A  similar  result  follows  from 
an  approximate  model  based  on  Poiseuille  flow,  described  in  Appendix  B.  The  observed 
particle  motion  was  therefore  not  attributable  to  absorption  of  optical  momentum  by  the 
liquid,  but  was  instead  consistent  with  thermophoresis.  Independent  simulations  of  natural 
(buoyant)  convection  were  also  performed  to  show  that  its  effect  was  negligible. 

3.6 


</> 

,sl 


0 

Fig.  7.  Simulated  z-velocity  of  the  liquid  in  sample  cell,  resulting  from  the  transfer  of  optical 
momentum  to  the  liquid  due  to  optical  absorption.  Velocities  shown  are  for  the  y  =  0  plane. 
The  black  square  marks  the  liquid-glass  boundary,  which  had  sides  50  pm  long. 


4.  Discussion 


At  this  point  a  comparison  between  the  thermal  diffusivity  DT  used  in  the  preceding 
calculation  with  that  of  various  theoretical  models  is  in  order.  For  this  purpose,  we  consider 
models  proposed  by  McNab  and  Meisen  [4],  Schimpf  and  Semenov  [3,6],  Ruckenstein  [8,9], 
Duhr  and  Braun  [10],  and  Brenner  [3],  for  which  the  thermal  diffusivities,  denoted  DTM,  DTS , 
Djri  Dtd  and  DTB ,  respectively,  are  given  by 

Dm= - r  a  i  3/7 - (12) 
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Table  2.  Material  Properties  Used  for  Model  Comparison 


Symbol 

Value 

Description 

P 

1033  kg/m" 

Density  of  liquid  [20] 

Too 

293  K 

Ambient  temperature 

kp 

1.35  W/m  K 

Thermal  conductivity  of  particle  [20] 

p 

7.0x10"*  K_1 

Volumetric  thermal  expansion  coefficient  of  liquid  [20] 

Cp 

2470  J/kg-K 

Isobaric  specific  heat  of  liquid  [20] 

n 

0.25  nm 

Molecular  radius  of  liquid 

AP 

6.5xlCT20  J 

Hamaker  constant  of  particle  [33] 

At 

5.6xlCT20  J 

Hamaker  constant  of  liquid  [34] 

£r 

29.3 

Relative  permittivity  of  liquid  [20] 

/pt 

-0.181  K 

[20] 

Idh 

1.1  pm 

Debye  length  of  liquid  (derived  from  [36,37]) 

Dtm 

5300(jm2/s-K 

Thermal  diffusivity  from  McNab-Meisen  model 

Dps 

0.012  |jm2/s-K 

Thermal  diffusivity  from  Schimpf-Semenov  model 

Dtb 

13  |jm2/s  K 

Thermal  diffusivity  from  Brenner  model 

Dtb 


k,P 


Pcp[1  +  {kp/2ki)_ 


(16) 


In  these  expressions,  p,  ft,  cp,  sR ,  rt  and  XDH  represent  the  liquid’s  density,  volumetric 
thermal  expansion  coefficient,  isobaric  specific  heat,  relative  permittivity,  molecular  radius, 
and  Debye  length,  respectively,  and  kp  is  the  particle’s  thermal  conductivity.  The  Hamaker 
constants  for  the  particle  and  liquid  are  denoted  as  Ap  and  A/.  Table  2  lists  the  relevant 
quantities  for  propylene  glycol  and  silica,  as  well  as  the  computed  thermal  diffusivities.  As  the 
interfacial  tension  of  the  propylene  glycol-silica  interface  y,  and  the  interaction  layer  thickness 
£  were  not  known  [1],  DTR  was  excluded  from  Table  2,  but  is  discussed  below.  As  the 
effective  surface  charge  per  unit  area  of  the  microsphere  oeff  was  also  unknown,  Dtd  is  also 
discussed  independently  below. 


4.1  The  McNab-Meisen  model 


The  McNab-Meisen  model,  which  is  simply  a  scaled  version  of  a  model  for  thermophoresis  in 
gases,  has  been  included  here  because  it  was  shown  to  accurately  model  the  thermophoretic 
velocity  of  individual  0.8- 1.0  pm  diameter  polystyrene  microspheres  in  water  [4].  As  shown 
in  Table  2,  the  predicted  value  for  DTM  exceeds  the  experimentally-fit  value  of22  pm2/s-K  by 
over  two  orders  of  magnitude,  so  this  model  could  not  be  regarded  as  consistent  with  the 
measured  data. 


4.2  The  Schimpf- Semenov  dipole-dipole  model 

The  Schimpf-Semenov  model  is  based  on  dipole-dipole  interactions  between  the  liquid  and 
particle,  and  was  originally  developed  for  particles  of  similar  size  to  the  liquid  molecules 
themselves  [3].  However,  it  was  extended  to  larger  particles  by  Brenner  [3]  in  his  derivation 
of  the  modified  form  of  Eq.  (13).  In  Table  2  it  is  shown  that  the  predicted  value  DTS  was  over 
two  orders  of  magnitude  lower  than  the  experimentally-fit  value  of  22  pm2/s-K.  Thus,  the 
Schimpf-Semenov  model  was  inconsistent  with  the  thermophoretic  motion  observed  in  the 
experiment. 

4.3  The  Ruckenstein  model 

The  Ruckenstein  model  is  based  on  the  presumption  that  temperature -induced  variations  in 
the  interfacial  tension  at  a  liquid-solid  interface  drive  thermophoresis  [8,9].  Although  data  for 
the  interfacial  tension  at  a  propylene  glycol-silica  interface  was  not  known  to  be  available  in 
the  literature,  the  fact  that  the  contact  angle  between  propylene  glycol  and  silica  is  small 
allowed  the  interfacial  tension  to  be  approximated  as  the  negative  of  propylene  glycol’s 
surface  tension.  Using  data  from  [35],  this  led  to  the  result  dy/dT  »7xl0"5  N/m-K  for 
propylene  glycol.  In  order  to  match  DTR  to  the  experimentally-fit  value  22  pm2/s-K,  the 
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interaction  layer  thickness  £  then  needed  to  be  26  nm.  By  way  of  comparison,  if  the  same 
procedure  were  performed  for  water  and  silica,  using  the  measured  thermal  diffusivity  22 
pm2/s-K  for  silica  in  pure,  deionized  water  [14],  and  the  interfacial  tension  dy/dT  «1.6xl0-4 
N/m-K  [20],  the  interaction  layer  thickness  would  instead  have  been  0.18  nm.  Given  the 
implausibility  that  the  interaction  layer  thickness  for  interfacial  tension  should  vary  so 
substantially  between  the  two  similar  liquids,  the  Ruckenstein  model  was  not  regarded  as  a 
viable  explanation  for  the  observed  thermophoretic  motion. 

4.4  The  Duhr-Braun  ionic -shielding  model 

The  Duhr-Braun  model  ionic -shielding  model  is  based  on  the  concept  that  under  conditions  of 
local  thermodynamic  equilibrium  and  constant  pressure,  steady-state  thermo -diffusion  will  be 
directly  related  to  the  entropy  of  the  particle-solvent  system  [10,15].  For  liquids  with 
sufficiently  weak  ionic  concentration  (such  as  pure  propylene  glycol),  this  leads  to  a  thermal 
diffusivity  that  is  proportional  to  the  Debye  length.  However,  in  order  to  satisfy  conditions  of 
local  thermodynamic  equilibrium,  it  is  necessary  that  the  temperature  gradient  be  limited  to 
the  range  [10] 


vr  < 


67rrjR2DT 


(17) 


or  equivalently,  through  Eq.  (10),  that  the  thermophoretic  velocity  be  limited  to 


\VT 


< 


kBT  ^ 

67irjR 2  ’ 


(18) 


where  kB  is  Boltzmann’s  constant.  For  conditions  in  the  experiment,  these  limits  were  |  VT\  < 
69  K/m  and  \vT\  <  0.0015  pm/s,  respectively,  or  nearly  two  orders  of  magnitude  below  those 
measured/simulated.  Thus,  a  treatment  of  thermophoresis  based  on  a  local  thermodynamic 
equilibrium  was  not  appropriate  given  the  conditions  of  the  experiment. 

4.5  The  Brenner  model 

The  Brenner  model,  in  its  most  basic  form  [11],  is  based  on  a  proposed  (yet  experimentally 
unconfirmed)  thermal  creep  boundary  condition  at  the  liquid-particle  interface,  which  directly 
leads  to  particle  thermophoresis.  The  expression  in  Eq.  (16)  also  embodies  a  novel 
interpretation  of  the  nature  of  viscous  stress  in  non-isothermal  fluids  [3].  Table  2  shows  that 
the  predicted  value  DTB  was  within  a  factor  of  1.8  of  the  experimentally-fit  value  22pm2/s-K. 
This  small  discrepancy  could  be  attributed  to  the  non-uniformity  of  VT0  in  the  experiment,  or 
inaccuracy  of  the  absorption  coefficient  2aj  and/or  heat  transfer  coefficient  h  used  in  the 
thermal  model.  It  also  may  be  due  to  the  fact  that  Eq.  (16)  did  not  account  for  “generation  of 
volume”  in  the  experiment  due  to  optical  absorption  (for  a  discussion  of  this  effect  see  [38]). 
Nonetheless,  the  relatively  good  agreement  demonstrates  a  sound  theoretical  basis  for  the 
thermophoretic  calculations  presented  in  Section  3.  Additionally,  these  results  provide  support 
for  the  Brenner  model  of  thermophoresis,  which  has  not  been  rigorously  tested  due  to  the 
scarcity  of  experimental  data  for  thermophoretic  motion  of  individual  particles  in  a  liquid 
[3,5]. 

5.  Extension  to  small  particles  in  an  optical  trap 

In  the  context  of  optical  trapping  and  manipulation  of  small  particles  in  liquids,  these  results 
are  significant  in  that  they  demonstrate  that  laser-induced  thermophoresis  can  under  certain 
circumstances  dominate  the  effects  of  radiation  pressure  in  a  conventional  optical  trap. 
Furthermore,  they  support  a  theoretical  model  that  can  be  used  to  quantify  the  impact  of 
thermophoresis  for  individual  particles.  In  order  to  expand  upon  these  points,  we  consider  an 
effective  thermophoretic  force  acting  on  a  spherical  particle  in  a  temperature  gradient,  defined 
in  analogy  to  Eq.  (4),  as 
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FTeq=67niRvT.  (19) 

This  quantity  is  referred  to  as  an  “effective”  force  because  it  does  not  necessarily  represent 
the  physics  behind  thermophoresis.  It  does,  however,  accurately  describe  the  thermophoretic 
contribution  to  the  net  force  on  a  particle,  under  the  common  assumption  that  the  net  force  is 
related  to  the  particle  velocity  according  to  Stokes’  law.  By  combining  Eqs.  (10)  and  (19),  the 
effective  thermophoretic  force  can  be  related  directly  to  the  temperature  gradient, 


F Teq  =  -67rrjRDTVT0.  (20) 

This  quantity  is  predicted  from  Eq.  (16)  to  scale  linearly  with  the  particle  radius  R.  Such 
scaling  is  also  supported  by  a  body  of  work  on  thermo -diffusion,  which  indicates  that  DT  is  for 
the  most  part  independent  of  R  in  the  quasi-hydrodynamic  limit  [1]  (with  the  notable 
exception  of  strongly-electrolytic  solutions  [10]).  In  contrast,  radiation  pressure  produces  a 
force  which  in  the  small  particle  limit  is  given  by  [17] 
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(  2  2  \ 
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(21) 


and  which  scales  at  least  as  R3.  Thus,  with  decreasing  particle  size,  thermophoresis  can 
eventually  dominate. 

In  the  simulations  presented  in  Section  3,  the  intensity  gradient  in  the  z-direction  had  peak 
value  0.105  W/pm3.  At  the  same  position  the  temperature  gradient  was  0.041  K/pm.  Using 
these  values  as  an  example,  along  with  material  properties  already  provided,  the  magnitudes 
of  both  the  effective  thermophoretic  force  and  the  maximum  force  due  to  radiation  pressure 
were  calculated  for  propylene  glycol  and  silica,  from  Eqs.  (20)  and  (21),  and  DT  =  22 
pm2/s-K,  and  plotted  versus  particle  radius  in  Fig.  8(a).  This  figure  also  plots  the  force  due  to 
radiation  pressure  for  a  polystyrene  particle  with  refractive  index  1.58.  Also  included  are  the 
effective  thermophoretic  forces  for  power  absorption  coefficient  of  19  m-1,  which  closely 
matches  that  of  both  propylene  glycol  and  water  at  1064  nm,  and  the  power  absorption 
coefficient  4900  m”1,  matching  that  of  water  at  1480  nm.  Note  that  the  latter  two  curves 
provide  a  comparison  of  the  relative  impact  of  thermophoresis  and  radiation  pressure. 
However,  in  practice,  temperature  constraints  may  limit  the  optical  power  that  can  be  used,  in 
which  case  both  forces  would  decrease  linearly  with  P. 

Figure  8(a)  predicts  that  for  silica  particles  suspended  in  propylene  glycol  (under  stated 
conditions),  thermophoresis  will  dominate  the  axial  radiation  pressure  throughout  the  nano¬ 
scale  (R  <  100  nm).  For  a  Gaussian  beam,  the  effective  thermophoretic  force  is  directed  away 
from  the  focus,  such  that  thermophoresis  serves  to  inhibit  trapping  in  3D.  Thus,  the  ability  (or 
inability)  of  a  Gaussian  beam  to  trap  small  particles  depends  not  only  on  random  thermal 
motion,  governed  by  kBT,  but  also  on  the  effective  thermophoretic  force. 
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Fig.  8.  Predicted  values  of  the  peak  radiation  pressure  force  and  the  effective  thermophoretic 
force,  along  the  beam  axis,  for  an  individual  silica  nano-sphere  in  pure  (a)  propylene  glycol  and 
(b)  deionized  water.  For  the  range  of  particle  radii  shown,  the  thermophoretic  force  dominates 
in  propylene  glycol.  For  deionized  water  the  thermophoretic  force  is  substantially  reduced,  but 
still  dominates  for  smaller  particle  radii. 

For  the  purpose  of  comparison,  Fig.  8(b)  plots  similar  curves  for  a  silica  nano-sphere 
immersed  in  pure,  deionized  water,  calculated  using  published  values  for  the  thermal 
diffusivity,  22  pm2/s-K  [14],  refractive  index,  1.33  [20],  and  dynamic  viscosity,  0.00089 
kg/m-s  [20].  For  simplicity,  the  intensity  and  temperature  gradients  used  in  these  calculations 
were  identical  to  those  used  for  propylene  glycol  (0.105  W/pm3  and  0.041  K/pm, 
respectively),  although  both  would  have  been  slightly  different  in  water.  As  shown  in  Fig.  8, 
the  effective  thermophoretic  force  is  predicted  to  be  much  weaker  in  deionized  water  than  in 
propylene  glycol.  This  is  a  consequence  of  the  fact  that  the  thermal  diffusivities  of  the  two 
liquids,  both  22  pm2/s-K,  are  remarkably  similar.  According  to  Eq.  (20),  for  an  identical 
particle  radius  and  temperature  gradient,  the  effective  thermophoretic  force  should  thus  be 
roughly  sixty  times  larger  in  propylene  glycol  than  in  deionized  water  (due  to  the  relatively 
low  viscosity  of  water).  That  is,  the  thermophoretic  force  scales  with  viscosity,  since  the 
thermal  diffusivity,  apparently,  does  not.  Such  a  concept  may  seem  counterintuitive,  but 
would  appear  to  be  consistent  with  Brenner’s  theories  on  the  nature  of  viscous  stress  in  a  non- 
isothermal  fluid,  and  thermophoresis  in  general  [3,11,38]. 

In  its  practical  application,  it  should  be  noted  from  Fig.  8(b)  that  the  thermophoretic  force 
on  a  silica  sphere  in  pure,  deionized  water  can,  but  will  not  necessarily,  dominate  the  force 
due  to  radiation  pressure.  For  higher  viscosity  liquids,  on  the  other  hand,  the  relative  strength 
of  the  effective  thermophoretic  force  offers  potential  advantages  for  the  control  of  individual 
particles  on  the  nano-scale. 

6.  Conclusion 

This  paper  has  presented  a  detailed  investigation  of  the  motion  of  individual  micro -particles  in 
a  moderately-viscous  liquid  in  direct  response  to  a  local,  laser-induced  temperature  gradient. 
The  elevated  viscosity  of  the  liquid  provided  for  substantial  differences  in  the  behavior 
predicted  by  various  models  of  thermophoresis,  which  in  turn  allowed  measured  data  to  be 
most  appropriately  matched  to  a  model  proposed  by  Brenner.  Based  on  this  model,  we  have 
quantitatively  predicted  when  thermophoresis  will  strongly  inhibit  the  ability  of  radiation 
pressure  to  trap  nano-scale  particles.  Furthermore,  the  model  predicts  that  the  thermophoretic 
force  will  scale  linearly  with  the  liquid  viscosity,  such  that  choice  of  liquid  plays  a  key  role  in 
the  relative  strength  of  the  thermophoretic  and  radiation  forces. 

Appendix  A 

We  calculate  the  axial  trapping  force  on  the  particle  using  the  method  outlined  by  Nieminen  et 
al  in  [24].  Briefly,  the  incident  beam  is  written  as  a  sum  of  vector  spherical  wavefunctions 
(VSWFs)  where  the  coefficients  of  the  VSWFs  that  correspond  to  the  incident  Gaussian  beam 
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are  calculated  using  a  point-matching  procedure  in  the  far  field  [25].  First,  the  far  field  limit 
for  a  paraxial  y-polarized  Gaussian  beam  as  a  function  of  the  beam  waist  is  used  to  estimate 
the  field  at  a  defined  set  of  points.  The  VSWFs  are  then  calculated  at  each  of  those  points  and 
the  resulting  overdetermined  system  of  equations  is  solved  in  a  least-squares  sense  for  the 
coefficients  of  the  VSWF  basis. 

Once  the  VSWF  coefficients  are  determined,  forces  on  the  particle  are  calculated  using  a 
T-matrix  approach  [26],  where  the  expansion  coefficients  of  the  scattered  wave,  P,  are  written 
as  a  function  of  the  incident  beam  coefficients,  A,  and  the  T-matrix  as  P=TA.  For  a  spherical 
particle,  the  T-matrix  is  diagonal  with  elements  given  by  the  usual  Mie  coefficients.  Given  the 
incident  and  scattered  coefficients,  the  force  and  torque  on  the  particle  can  be  computed  [24]. 

We  note  that  the  incident  coefficients  corresponding  to  a  single  origin  and  orientation  need 
only  be  calculated  once  using  the  point-matching  procedure.  Once  these  have  been 
determined,  the  coefficients  corresponding  to  other  spatial  locations  in  the  same  trapping  field 
can  be  determined  using  coordinate  transformations  and  rotations  of  the  VSWFs  [24].  We  first 
determine  the  coefficients  for  a  VSWF  basis  centered  at  the  beam  waist  and  then  predict  the 
forces  corresponding  to  particle  locations  at  various  points  along  the  axis  of  the  beam. 

Appendix  B 

From  Eq.  (11),  the  optical  momentum  absorbed  by  the  liquid  produces  a  force  per  unit  length 
across  the  cross-section  of  the  optical  beam  of 


dz 


=  J  {  fal,sdxdy  = 


2aPntz 


(Bl) 


In  comparison,  the  laminar  flow  of  liquid  through  a  circular  channel  of  radius  achanneh 
assuming  no-slip  boundary  conditions,  is  related  to  the  force  per  unit  length  across  the  channel 
[39] 


~^  =  ^{^achonnel)z=^flVmm,  (B2) 

dz  dz  ' 

where  p  is  the  fluid  pressure  and  vmax  the  peak  velocity  of  the  liquid.  Equating  Eqs.  (Bl)  and 
(B2)  leads  to  a  rough  estimate  of  the  fluid  flow  resulting  from  the  absorption  of  optical 
momentum  by  the  liquid 


aPn^z, 

27TTJC 


(B3) 


When  related  to  experimental  conditions,  this  predicts  a  peak  velocity  of  0.0007  pm/s,  which 
was  within  a  factor  of  5  of  that  determined  by  the  more  rigorous  finite  element  simulations. 
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